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ABSTRACT 

The observed spectra of 9 pulsars for which multiwavelength data are available from 
radio to X- or 7-ray bands (Crab, Vela, Geminga, B0656+14, B1055-52, B1509-58, 
B1706-44, B1929+10, and B1951+32) are compared with the spectrum of the radia- 
tion generated by an extended polarization current whose distribution pattern rotates 
faster than light in vacuo. It is shown that by inferring the values of two free pa- 
rameters from observational data (values that are consistent with those of plasma 
frequency and electron cyclotron frequency in a conventional pulsar magnetosphere) , 
and by adjusting the spectral indices of the power laws describing the source spectrum 
in various frequency bands, one can account quantitatively for the entire spectrum of 
each pulsar in terms of a single emission process. This emission process (a generaliza- 
tion of the synchrotron-Cerenkov process to a volume-distributed source in vacuum) 
gives rise to an oscillatory radiation spectrum. Thus, the bell-shaped peaks of pulsar 
spectra in the ultraviolet or X-ray bands (the features that are normally interpreted as 
manifestations of thermal radiation) appear in the present model as higher-frequency 
maxima of the same oscillations that constitute the emission bands observed in the 
radio spectrum of the Crab pulsar. Likewise, the sudden steepening of the gradient 
of the spectrum by —1, which occurs around 10 18 — 10 21 Hz, appears as a universal 
feature of the pulsar emission: a feature that reflects the transit of the position of the 
observer across the frequency-dependent Rayleigh distance. Inferred values of the free 
parameters of the present model suggest, moreover, that the lower the rotation fre- 
quency of a pulsar, the more weighted towards higher frequencies will be its observed 
spectral intensity. 

Key words: radiation mechanisms: non-thermal — pulsars: individual: Crab pulsar, 
Vela pulsar, Geminga pulsar, PSR B0656+14, PSR B1055-52, PSR B1509-58, PSR 
B1706-44, PSR B1929+10, PSR B1951+32. 



1 INTRODUCTION 

Pulsars are unique in emitting a radiation whose spectrum 
extends from radio waves to gamma rays with essentially 
the same characteristics ( |Lyne &; Graham-Smith|20 06 ) . The 



fact that in many pulsars the pulses emitted over widely 
separated frequency bands (e.g. in the radio and X-ray 
bands) are closely aligned in phase and have correlated pro- 



files ( |Moffett fe Hankins||l996| |Kuiper et al.||2001[ |2002 
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Cusumano et al.| [20Q3; Lo mmen et al.||2007| indicates that 
both the mechanism by which such pulses are generated and 
the magnetospheric site from which they originate must be 
the same for all frequencies. This notwithstanding, the ex- 
tant models of pulsar radiation invoke a diverse set of unre- 
lated emission mechanisms and emission sites to account for 
the observational data on multiwavelength characteristics of 
this radiatio n (|Ruderman &; Sutherland! 1975| |Caraveo et al.| 
2004| |Zavlin||2007| |Harding et al.||2008|). 

It was shown in |Ardavan et al.| ( |2008a ) that the rigid ro- 
tation of the overall distribution pattern of the pulsar emis- 
sion can only arise from an emitting current whose distri- 
bution pattern likewise rotates rigidly with the same an- 
gular frequency (see also Section 2). The fact that this re- 
quires the pattern of distribution of the emittting current 
to move with a linear speed exceeding the speed of light in 
vacuo outside the pulsar's light cylinder is not incompatible 
with special relativity: the charge separation resulting from 
the coordinated motion of an aggregate of charged parti- 
cles gives rise to a polarization current whose distribution 
pattern can have a superlumi nal motion (|Ginzburg||1972| 
|Bolotovskii fe Ginzbur"g||1972| |Bolotovskii fe Bykov||1990| ). 
Moreover, Maxwell's equations show that a polarization cur- 
rent contributes towards the radiation field in just the same 
way as a current of free charges. In fact, such superluminal 
sources of radiation have been experimentally realized by 
several groups ([Bessarab et al.||2004| |Singleton et al.| |2004 ; 



Ardavan et al.|2004b| |Bolotovskii fe Serov|2005| ). 

Here we show that the spectrum of the emission from 
a polarization current whose distribution pattern rotates 
faster than light in vacuo fits the spectra of 9 pulsars, 
for which multiwavelength data are available, over the en- 
tire range of their output from radio waves to gamma 
rays. The model we adopt for the source is a generic one: 
it has an azimuthally-fluctuating distribution pattern that 
both rotates and oscillates [as realized in recent experi- 
ments ( |Singleton et al.|[2004| |Ardavan et al.||2004b| ], and 
it has sharp spatial fluctuations in the direction parallel 
to the rotation axis with power-law spectra in various fre- 
quency bands [as demanded by recent numerical models of 
the pulsar magnetosphere that predict the formation of cur- 



rent sheets (Spitkovsky 2006)]. We infer the values of the 
two adjustable parameters, characterizing the frequencies of 
temporal and azimuthal fluctuations of the source, from the 
observational data, and choose the spectral indices of the 
power laws describing the spectral distribution of the varia- 
tions of the source parallel to the rotation axis by optimizing 
the goodness of fit to the data. We shall see that the the- 
oretical spectra thus derived from a single emission process 
(a process that takes place in a localized region of the mag- 
netosphere just outside its light cylinder) can quantitatively 
account for the data on all nine of these extensively observed 
pulsars over 15 to 18 orders of magnitude of frequency. 

The emission mechanism in question may be regarded 
as a generalization of synchrotron- Cerenkov process to a 
volume-distributed source in vacuo. Its spectral distribution 
is described by the square of a Bessel function whose argu- 
ment exceeds its order, and so has an oscillating amplitude 
with an algebraic (rather than exponential) rate of decay 
with frequency ( Ardavan et al.|2 003). Giving rise to a radi- 
ation whose spectrum is oscillatory, it is a mechanism that 
naturally accounts for ( Ardavan et al.|2008a| the occurrence 



of the observed emission bands in the dynamic spectrum of 
the Crab pulsar ( |Hankins fe Eilek||2007| ). 

The unified explanation this mechanism offers for the 
apparently unrelated features of pulsar spectra in different 
spectral bands sharply contrasts with the variety of models, 
proposed in the published literature, which could at best ex- 



plain each of these features separately ( Ruderman & Suther- 



land|T975| [Caraveo et al.|2004| [Zavlin||2007| |Harding et a! 



2008). The bell-shaped peak of the spectrum in the ultravio- 



let or X-ray bands that is normally attributed to black-body 
radiation from hot spots in the pulsar magnetosphere ( Car- 
|aveo et al.|2004| ), for instance, emerges as a higher-frequency 
maximum of the same proportionately-spaced oscillations of 



the spectrum that Hankins & Eilek (2007) have observed in 



the radio band. Moreover, the decrease in the spectral index 
of the radiation by —1 that is universally encountered in 
pulsar spectra around 10 18 — 10 21 Hz turns out to arise, not 
from any changes in the spectral distribution of the source 
of the radiation, but from a transition though the Rayleigh 
distance (see Section [3]). 

This paper is organized as follows. In Section |2.1[ we 
point out why the observational data on pulsars demand 
that the source of pulsar radiation should have a superlu- 
minally rotating distribution pattern. In Section 2.2, we de- 
scribe the field due to a constituent volume element of such 
a source, i.e. the Lienard-Wiechert field generated by the su- 
perluminal counterpart of synchrotron process ( | Schmidt et| 
aLl|2007| jArdavan et aLl|2004a| ) . In Section [53] we consider 



the field arising from the entire volume of the source. We 
point out how the most efficient parts of a pulsar magneto- 
sphere for the formation of signals detectable at large dis- 
tances are the thin filaments within the superluminally ro- 
tating part of its current distribution pattern that approach 
the observer with the speed of light and zero acceleration at 
the retarded time ( Ardavan et al.|[2007| . Because the wave 
fronts from each such filament form a caustic on which the 
contributions that are made towards the field over a finite 
interval of emission time are received during a considerably 
shorter interval of observation time, we will see that the 
overall radiation beam generated by a volume source con- 
sists of an incoherent superposition of coherent, narrowing 
subbeams whose intensities diminish as i?p _1 , i nstead of 
Rp~ 2 , with the distance Rp from their source (Ardavan 



et al.|[2004a ). Consequently, the radiation detected at large 



distances from the pulsar will always be dominated by these 
subbeams. In Section [3] we fit the multiwavelength obser- 
vational data on each of the 9 pulsars listed above with the 
frequency spectrum predicted for a rotating superluminal 
source; our earlier derivation of the spectrum (Ardavan et| 



al. 2003) is outlined from an alternative point of view in 



Appendix A. Implications of the inferred values of the pa- 
rameters of the fitted spectra will be discussed in Section [4] 
and a summary will be given in Section [5] 



2 SUPERLUMINAL MODEL OF PULSARS 

2.1 Observational constraints on the motion of 
the distribution pattern of the source 

The rigid rotation of the overall distribution pattern of the 
pulsar emission reflects a radiation field E whose cylindrical 
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components depend on the cylindrical coordinates (r, ip, z) 
and time t as 



z (r,<fi,z;t) = E r ^ z (r,ip -ut,z,t), 



(1) 



where oo is the rotation frequency of the pulsar. Such a field 
can only arise from an electric current whose density j like- 
wise depends on the azimuthal angle ip in the combination 
(p — u)t only: 



j r ,<p,z (r, <p, z\ t) = j r ,<p,z (r, ip - out, z, t) 



(2) 



[see Appendixes A and B of Ardavan et al. (2008a)]. This 
property of the emitting current follows not only from 
the observational data, but also from the numerical mod- 
els of the magnetospheric structure of an oblique rota- 
tor; it is found that any time-dependent structures in such 
models rapidly approach a steady state in the corotating 
frame flSpitkovsky1|2006 ). 

Unless there is no plasma outside the light cylinder, 
therefore, the emitting polarization current distribution in 
the magnetosphere of a pulsar should have a superlumi- 
nally rotating pattern for r > c/co (where r is the radial 
distance from the axis of rotation and c is the speed of 
light in vacuo). Such a source is not inconsistent with spe- 
cial relativity, as the charged particles that make up the 
polarization current distribution need only move relatively 
slowly ( |Ginzburg|1972||Bolotovskii & Ginzburg|1972| ). It has 
been experimentally verified that such moving charged pat- 
terns act as sources of radiation in precisely the same way as 



any other moving sources of electromagnetic fields ( Bessarab 



et al.|2004| | Ardavan et al.|2004b| [Singleton et al.|2004||Bolo- 
tovskii & Serov||2005| 



Many distinctive features of the emission from a super- 
luminal source are illustrated by the radiation from the fol- 
lowing generic polarization current whose distribution pat- 
tern rotates and oscillates at the same time: j = dP/dt for 
which 



Pr,<p,z(r,<p,z,t) 



s r ,<p,z(r, z) cos(m(f) cos(Ot), 

— 7T < (f < 7T, (3) 

and op = op — out. Here, P r ,<p,z are the components of the 
polarization P in a cylindrical coordinate system based on 
the axis of rotation, s(r, z) is an arbitrary vector function 
with a finite support in r > c/w, m is a positive integer, 
and O is an angular frequency whose value differs from an 
integral multiple of the rotation frequency oo [for the signifi- 
cance of this incommensurablity requirement, see Appendix 



A and |Ardavan et al. (2003)]. For a fixed value of t, the 
azimuthal dependence of the polarization ^ along each cir- 
cle of radius r within the source is the same as that of a 
sinusoidal wave train with the wavelength 2nr /m whose m 
cycles fit around the circumference of the circle smoothly. 
As time elapses, this wave train both propagates around 
each circle with the velocity roo and oscillates in its ampli- 
tude with the frequency Q. This is a generic source: one can 
construct any distribution with a uniformly rotating pat- 
tern, P rj(p7Z (r,(p, z), by the superposition over m of terms 
of the form s r ,cp, z (r, z, m) cos(ra<£). It also corresponds to 
laboratory-based sources that have been used in experimen- 
tal demonstrations of some of the phenomena described be- 
low ( |Ardavan et al.|2004b||Singleton et al.|2004 ). 

The results reported here are derived from the retarded 
solution of Maxwell's equations for the above current dis- 




(a) 




Figure 1. (a) Envelope of the spherical wave fronts emanating 
from a superluminally moving source element (S) in uniform cir- 
cular motion. The heavy curves show the cross section of the 
envelope with the plane of the orbit of the source. The larger of 
the two dotted circles designates the orbit and the smaller the 
light cylinder, (b), (c) and (d) Space-time (ct versus x) diagrams 
showing the intersection of the trajectory of the source point S 
with the past light cone of the observation point P when P lies 
outside (b), inside (c), and on the cusp of (d) the envelope of wave 
fronts. 



tribution [see I Ardavan et al.| j2008a| |2007[ |2004a[ |2003 ) 
Schmidt et al.f fl2007|); |Ardavan| |p998f]. 



2.2 The field generated by a single volume 
element of the source 

A superluminal source is necessarily 



distributed (Bolotovskii & Ginzburg 1972). 



volume- 
Howe ver, 

its field can be built up from the superposition of the 
fields of its moving constituent volume elements which are 
essentially point-like. Figure [TJa) shows that the waves 
generated by a constituent volume element of a rotating 
superluminal source possess a cusped envelope and that, 
inside the envelope, three wave fronts pass through any 
given observation point simultaneously. This reflects the 
fact that the field inside the envelope receives simultaneous 
contributions from three distinct values of the retarded 
time [see Fig. [TJc)]. On the cusp of the envelope, where the 
space-time trajectory of the source is tangent to the past 
light cone of the observer [Fig.[TJd)], all three contributions 



toward the value of the field coalesce (|Ardavan| [1998 
Ardavan et al.|2003| |2004a| ) 



On this cusp (caustic), the source approaches the ob- 
server with the speed of light and zero acceleration at the re- 
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Figure 2. Three dimensional views of the envelope (a) and its 
cusp (b). 



tarded time, i.e. dR(t)/dt = —c and d 2 R(t)/dt 2 = 0, where 
R(t) = |x(t) — xp| is the distance between the source point 
x(t) and the observation point xp. As a result, the interval 
of emission time for the signal carried by the cusp is much 
longer than the interval of its reception time ( Ard avan et al.| 



2003). This essentially instantaneous reception of contribu- 



tions from an extended period of emission time represents 
focusing of the radiation in the time domain: a unique effect 



that has already been demonstrated experimentally (Arda- 
van et al.|2004b| | Singlet on et al.|2004) 



A three-dimensional view of the envelope of wave fronts 
and its cusp is shown in Fig. [2] The two sheets of the en- 
velope, and the cusp along which these two sheets meet 
tangent ially, spiral outward into the far zone. In the far 
zone, the cusp lies on the double cone Op = arcsin[c/(ra;)], 
Op = 7T — arcsin[c/(rcj)], where (Rp,6p, ipp) denote the 
spherical polar coordinates of the observation point P. Thus, 
a stationary observer in the polar interval arcsin[c/ (rco)] ^ 
Op ^ 7r — arcsin[c/(ro;)] receives recurring pulses as the en- 
velope rotates past him/her ( | Ardavan et al.|2 007). 

Figure [3] shows the radiation field generated by the ro- 
tating source element S on a cone close to the cusp, just 
outside the envelope. Not only does the spiralling cusp em- 
body a recurring pulse, but the plane of polarization of the 
radiation swings across the pulse ( |Schmidt et al.|2007|, as in 



the radio emission received from pulsars ( |Lyne &; Graham- 
Smith|[2006 ) 



As a consequence of the multivaluedness of the retarded 
time, three images of the source are observeable inside the 
envelope at any given observation time (Fig. [4|. The waves 
that were emitted when the source was at the retarded posi- 
tions ii, I2 and I3 in Fig.[4]are all received simultaneously at 
the observation point P. These images are detected as dis- 
tinct components, or modes, of the emitted radiation ( |Bolo-| 
|tovskii fe Bykov|1990 ). 

Field strengths and polarization position angles of the 
three images (radiation modes) close to the cusp are shown 
in Fig. [5] Two modes dominate everywhere except in the 
middle of the pulse. Moreover, the position angles of two 
of the modes are 'orthogonal' and that of the third swings 
across the pulse bridging the other two. The constructive in- 
terference of the emitted waves on the envelope (where two 
of the contributing retarded times coalesce) and on its cusp 
(where all three of the contributing retarded times coalesce) 
gives rise to the divergence of the field of a point-like source 
on these loci ( Ardavan|1998 ). Here we plot the spatial distri- 
bution of the field excising the narrow regions in which the 



3 a. 



-5 




-5 



0, 
x p co/c 



Figure 3. Polarization position angles and field strengths on the 
cone Op = 7r/12 outside the envelope of wave fronts for a source 
element with ruj = 2. Here, the field vectors are projected onto 
the plane (xp,yp) of the source's orbit. 




Figure 4. An observer P inside the envelope detects three images 
(ii, I2, 1 3) of the source S simultaneously. 



certain threshold (Schmidt 



magnitude of the field exceeds 
|et al.||2007| ). 



2.3 The field generated by the entire volume of 
the source 

2.3.1 Nondiffracting subbeams comprising the overall 

beam and the nonspherical decay of their amplitudes 

The dominant contribution towards the field of an extended 
source comes from a thin filamentary part of the source 
that approaches the observer, along the radiation direction, 
with the speed of light and zero acceleration at the retarded 
time ( Ardavan et al.||2007 ). For an observation point P in 
the far zone with the coordinates (Rp,0p, (fp), this filament 
is located at r = (c/00) esc Op, ip = ipp + 37r/2 and is essen- 
tially parallel to the rotation axis (Fig. (6|. The collection 
of cusps of the envelopes of wave fronts that emanate from 
various volume elements of the contributing filament form 
a subbeam whose polar width is nondiffracting: the linear 
dimension of this bundle of cusps in the direction paral- 
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Figure 6. Schematic illustration of the light cylinder, the filamen- 
tary part of the source that approaches the observation point with 
the speed of light and zero acceleration at the retarded time, the 
orbit r = c/(lj sin Op) of this filamentary source, and the subbeam 
formed by the bundle of cusps that emanate from the constituent 
volume elements of this filament. 
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Figure 5. (a) The relative strengths of the three radiation modes 
as observed near the cusp on a sphere of large radius. The to- 
tal field strength (black) and strengths of the underlying contri- 
butions from the three images of the source (green, red, blue) 
are shown for a source element with ruj = 1.1 and an observa- 
tion point that sweeps a small arc of the circle Rpuo/c = 10 10 , 
$P = 7r/2.7, crossing the envelope near the cusp, (b) The cor- 
responding position angles of the contributions from the three 
retarded times (green, red, blue) are shown relative to one an- 
other and to that of the total field (black); the position angles of 
the dominant contributions are shown with open triangles, and 
those of the weaker contributions with filled triangles. 



lei to the rotation axis remains the same at all distances 
from the source, so that the polar angle 50 p subtended by 
the subbeam decreases as Rp -1 with increasing Rp [Fig. M 



and |Ardavan et al.| ([2007])] . 

In that it consists of caustics and so is constantly dis- 
persed and reconstructed out of other waves, the subbeam 
in question radically differs from a conventional radiation 
beam [see Appendix D of Ardavan (1998)]. The narrow- 
ing of its polar width (&s Rp~ ) is accompanied by a more 
slowly diminishing intensity (an intensity that diminishes as 
Rp -1 instead of Rp~ 2 with distance), so that the flux of 
energy across its cross sectional area remains the same for 



all Rp (Ardavan et al. 2007). This slower rate of decay of 
the emission from a superluminally rotating source has been 



confirmed experimentally ( Ard avan et al.|[2004b| |Singleton 
et al.||2004| ) 



of light and zero acceleration) changes as the source rotates 
(see Fig. [6}. In the case of a turbulent plasma with a super- 
luminally rotating macroscopic distribution, therefore, the 
overall beam within which the narrow, nonspherically de- 
caying radiation is detectable would consist of an incoher- 
ent superposition of coherent, nondiflxacting subbeams with 
widely differing amplitudes and phases [similar to the train 
of giant pulses received from the Crab pulsar (Slowikowska 
et al.||2005D ]. 

The overall beam occupies a solid angle whose polar and 
azimuthal extents are independent of the distance Rp. It 
is detectable within the polar interval arccos[(r<o;/c) _1 ] ^ 
\0p — tt/2\ ^ arccos[(r>o;/c) _1 ], where [r<,r>] denotes the 
radial extent of the superluminal part of the source. The az- 
imuthal profile of this overall beam reflects the distribution 
of the source density around the cylinder r — c/(cjsin#p), 
from which the dominant contribution to the radiation 
arises ( |Ardavan et al.||2007| ). 

The fact that the observationally-inferred dimensions 
of the plasma structures responsible for the emission from 
pulsars are less than 1 metre in size ( |Hankins et~ld^|20(^ 
reflects, in the present context, the narrowing (as .Rp - and 
Rp~ 3 , respectively) of the radial and azimuthal dimensions 
of the filamentary part of the source that approaches the 
observer with the speed of light and zero acceleration at 
the retarded time (Ardavan et al. 2007). Not only do the 



nondiflxacting subbeams that emanate from such filaments 
account for the nanostructure, and so the brightness tem- 
perature, of the giant pulses, but the nonspherical decay of 
the intensity of such subbeams (as Rp -1 instead of Rp~ 2 ) 
explains why their energy densities at their source appear 
to exceed the energy densities of both the plasma and the 
magnetic field at the surface of a neutron star when esti- 



The contributing part of an extended source (the fila- 
ment that approaches the observation point with the speed 



mated on the basis of the inverse square law (Soglasnov et 
|aL|2004) . 



2.3.2 Frequency spectrum of the radiation by a volume 
source 

The spectrum of the radiation emitted by an extended 
source with a superluminally rotating distribution pattern 
is oscillatory with oscillations whose spacing increase with 
frequency [see Ardavan et al. (2008a, 2003} and equation Q 
below]. While the Bessel function describing synchrotron 
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Figure 7. The predicted oscillations of the spectrum of the emis- 
sion for uj/(2tt) = 30.3 Hz and fl/uu ~ 1.9 x 10 4 , shown in (a), have 
the same spacing as those of the emission bands in the observed 
spectrum of the Crab pulsar ( |Hankin~ &; Eilek 2007), shown in 
(b) [see |Ardavan et al.| ( |2008a| )]. 



radiation has an argument smaller than its order and so 
decays exponentially with increasing frequency, the Bessel 



function appearing in the present analysis (Ardavan et al. 



2003), whose argument exceeds its order, is an oscillatory 



function of frequency with an amplitude that decays only 
algebraically ( jAbramowitz &; St egun 1970). Figure [7] shows 
that the spacing of the emission bands in the spectrum of 
the Crab pulsar fit the predicted oscillations for an appro- 
priate choice of the single parameter Q/lj. The value of this 
parameter, thus implied by the data of Hankins & Eilek 
(20071, places the last peak of the oscillating spectrum at 
a frequency (~ Q 3 /to 2 ) that agrees with the position of the 
ultraviolet peak in the spectrum of the Crab pulsar. By in- 
ferring the remaining free parameter m in equation ^ from 
the observational data and by adjusting the spectral indices 
of the power laws describing the spectral decomposition of 
the source density s in various frequency bands [see Sec- 
tion (J3j)] , we are thus able to account for the continuum 
spectrum of the Crab pulsar over 16 orders of magnitude of 
frequency (Fig. [8). 

At radiation frequencies higher than ~ 10 18 Hz, the 



Rayleigh distance [see Section 3.5 and Brooker (2003)] be- 
comes larger than the distance of the Crab pulsar from the 
Earth. For an observer within the Rayleigh distance, the 
emission arises from a narrower radial extent of the source, 
and so has a higher degree of mean polarization, the higher 
the frequency. The degree of circular polarization of such a 
high-frequency emission decreases with increasing frequency 




logio v (Hz) 

Figure 8. The points show observational data (where available) 
on the spectrum of the Crab pulsar. The curves show the spectral 
distribution log(dP n /df2p), predicted by equations |8j and 
j9|, versus logn and log v for v = nuj /(2ir) ~ 30. 3n Hz and 
Q/cu c± 1.9 x 10 4 . In the model, the recovery of intensity at the 
ultraviolet peak (~ 10 15 Hz) is caused by resonant enhancement 
due to the azimuthal modulation frequency muo / '(2tt) ~3 x 10 13 
Hz. The steepening of the gradient of the spectrum by —1 at 
2.4 x 10 18 Hz reflects a transition through the Rayleigh distance. 



-1/3 



so that this emission is essentially 100 per cent 
linearly polarized at all phases, including off-pulse phases 



[see equation (A26) below]. Note that, in this model, the in- 
tensity and polarization of the off-pulse emission, too, reflect 
the distributions of density and orientation of the emitting 
current around the cylinder r — c / \uj sin p) within the pul- 
sar magnetosphere ( Ardavan et al.||2007 ). 

We will see in Appendix A that the reason the spec- 
trum of the emission from the source distribution Q is not 
limited, in the superluminal case, to just the frequencies 
Q ± moo is the following. On the one hand, each element up 
of a constituent ring of the source makes its contribution 
towards the observed field at the discrete set of retarded 
times Lot — ipp — <p -\- 3ir/2 + 2kir at which it approaches 
the observer with the speed of light and zero acceleration, 
i.e. periodically with the period 2tt /lo (k is an integer). On 
the other hand, temporal modulations of the density of this 
source occur on a time scale 2tt/Q that is incommensurable 
with the period 2tt/lo. Because the emission time of the ele- 
ment labelled by up is fixed (via cot = up P — up -\- 3tv /2 + 2kn) 
by its initial azimuthal position (op), the temporal [cos(Qt)] 
and spatial [cos(mup)] modulations of source (|3) effectively 
combine into a single variation [oc exp(iQup / 'lo) cos(ra<£)]. It is 
the Fourier decomposition of this variation with up, or equiv- 
alently with the retarded time t, that together with the in- 
commensurablity of the two periods (2ty/co and 27r/Q), and 
the finiteness of the domain of definition of up, results in a 
spectrum containing all frequencies [see Fig. 5 of Ardavan 
etaDd2003b]. 
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3 THEORETICAL FITS TO THE 

MULTIWAVELENGTH DATA ON PULSAR 
SPECTRA 

3.1 The equation used to create the fits 



Using the result derived in Appendix A [equation A30)], 
we can write the frequency dependence of the power that 
is radiated by the source distribution |3} per harmonic per 
unit solid angle as 



dP n 



5i(n) 2 Ai 2 



(!) 



!/ 3 fi 
0) 



+ &(»)'(§) 



2/3 



Ai' 

1/3 



./2 



Ai 



■(1) 
(!) 



" 3 Q 



x Ai' 
in which 

5i(n) = n 2/3 |K r |l^ 

&(n)=n 2/3 |K r |^ 
and 

S 3 (n) = 



(!) 



(4) 



Sr + Sep COS^ 



\sJ 2 + |s r | 2 cos 2 ( 



s z sin t 

\l/2 



,2\ 1/2 



2/3 |K,|K 



^ I s r cos 6 



(s^ cosOp — s z sin Op) — s r s* 



, 1/2 



Here, s r ,<p, z are the Fourier components of the source densi- 
ties s r ,ip,z\r=csc9 P with respect to z [see equation (A27)], 



S and the superscript star denote the imaginary part 
and the conjugate of a complex variable, respectively, and 
(Rp,9p,ipp) are the spherical polar coordinates of the ob- 
servation point P. The function K^ Q is defined by 



(-1) 



n+m 



+ 



(5) 



and the various frequencies involved (n, etc.) are tab- 
ulated for convenience in Table [I] Finally, the function K r 
is defined for various frequency ranges in Table [2] f< and 
r> denote the radial boundaries of the part of the source 
that contributes towards the radiation at P in units of the 
light-cylinder radius c/co. 



Table 1. Frequencies involved in the fitting of the pulsar spectra. 
For convenience in comparisons with observations, most param- 
eters are presented in frequency units (i.e. in Hz) rather than as 
angular frequencies (i.e. in radians per second). 



Frequency 


Definition 


10 


pulsar rotation frequency 




observation frequency 


mio 
2tt 


azimuthal modulation frequency 


n 

2tt 


temporal modulation frequency 




^ + m 


fl- 


n 

— — m 

CO 



Table 2. Limiting values of parameters and functions in various 
frequency ranges; the harmonic number n yields the observation 
frequency v = nuj/27T in units of the rotation frequency uj/2tt. 



Function 



Frequency range 



Limiting value 



/ 2 \ V 3 n 

V n ) to 



■«(! 



Ai 



2 \V3 n 



Ai(0) 



Ai' 



V3 n 



-«(er 



Ai' 



-»(S) £ 



Ai'(O) 



n <C /i+ 



(f>-f<) 2 



: r> - r< 



n > 



(f>-f<) 2 



(2zA P )'exp(- 



spectrum of the radiation described by equation Q is deter- 
mined by the ratio Q/uj. Observational data specify the value 
of this ratio for each of the 9 pulsars we consider quite accu- 



3.2 Fit procedure 

The pulsar spectra are fitted to the predicted spectrum Q 
by adjusting the modulation frequencies Q/2n and moo /2tv 
of the source (see Table [I]), and by specifying the required 



rately [see also |Ardavan et aL (2008a)]. However, the existing 
data are insufficient to determine the remaining parameters 
of the source [frequency muj of the azimuthal modulations of 
the current distribution Q and the frequency dependence of 
the Fourier components s rjlp , z of the amplitude of this cur- 
rates of change of the Fourier components s r ,<p,z of the source rent] uniquely. We will assume that the amplitudes s r ,<p,z [or 



density (or equivalently, and more practically, Si, S2 and 
S3) with frequency v — nuo/2n. Our objective in each case 
is to illustrate that the observed spectrum can be faithfully 
represented by that in equation Q within a wide margin 
of accuracy. We will model the frequency dependence of the 
factors that multiply the Airy functions in this equation by 
as simple a set of power laws as allowed by the available 
data. We will also extrapolate the spectra thus derived into 
the frequency intervals for which no data are available. 

The position of the peak of the last oscillation in the 



equivalently the functions Si(n), i = 1,2,3, that appear in 
equation Q] have power-law dependences of the form n ai 
on the harmonic number n with differing exponents oli in 
different frequency bands. The relative amplitudes, S2/S1 
and S3/S2, and the values of ai are then inferred from the 
data. 

We start plotting each figure with the inferred values of 
the above adjustable parameters in the radio band and con- 
tinue the resulting theoretical curves as far as the next patch 
of available data for which one or more of these parameters 
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have different values. In a few cases (Vela, Geminga and 
B0656+14), we will cover an interval for which no data are 
available by extrapolating the theoretical curves implied by 
the data on opposite sides of the empty interval until they 
meet. We are not predicting that the spectrum would neces- 
sarily follow the plotted curves in frequency intervals where 
no data are available. Not only is the polarization current 
distribution in the magnetosphere of a pulsar not expected 
to have such a simple space-time structure as that described 
by equation ([3]), but even in such a simplified description, 
the amplitudes s r ,^, z need not have the power-law forms we 
have approximated them by. 

Before embarking on a detailed description of the fits, 
we show, in the following three Subsections, how the be- 
haviour of the various components of equation Q (see Ta- 
ble [5]) leads to some "universal" features of the emission 
spectra of pulsars. Evidently, many of the gross phenomena 
seen in pulsar observations are quite general consequences 
of the fact that they arise from a rotating, superluminal 
source. It is only in modelling the detailed behaviour of the 
spectra that it is necessary to specify parameters Si , S2 and 
S3 (which represent intricacies of the 'weather' around the 
emitting region as opposed to the 'climate' in the pulsar 
magnetosphere) with any accuracy. 



3.3 Resonant enhancement 

According to equation ([5]), the value of K^ Q (and, hence, 
those of Si , S2 and S3) can undergo a sudden change as 
a result of the resonance between the radiation frequency 
a;n/27r and one of the source frequencies cj/j j ±/2tt. As sum- 
marized in Table [2] K<p is independent of n (~ 1) when 
n <C decays as n _1 when n ^> and is of the or- 
der of n when n ~ m ^> Q/oo [see equation^]; hence, its 
value sharply increases by several orders of magnitude for 
n rsj m ^> Q/uj. This is the primary cause of the recovery of 
the intensity at higher frequencies that will be seen in the 
fitted spectra below. 



3.4 The Airy function and its derivative 

The Anger functions appearing in the general expression 



for the radiated power [equation (66) of Ardavan et al. 
(2003)] reduce to Bessel functions for n ^> 1 and to Airy 
functions for n ^> Q/cj [see Appendix C of Ardavan et 



[ah] ( |2008a| ]. In most parts of the spectra modelled here, 
Q/lu <C n <C (Q/u) 3 and so Ai'/Ai is greater than unity (Ta- 
ble [2]). Consequently, when the second term of equation Q 
is employed in the fits below, S2 / Si is given values of the or- 
der of 10 6 . In other words, where invoked, the second term is 
actually comparable to the first term everywhere apart from 
the region past the last maximum. The dependence on n of 
band spacings and amplitudes of the two functions Ai and 
Ai' are sufficiently different for the inclusion of the second 
term to make a difference to the shape of the spectrum. 

The third term of equation Q declines with n at a 
similar rate to the second term in regions of the spectrum 
where n « (Q/00) 3 . However, it was not found necessary 
to invoke this term in the fits below, probably because S3 is 
much smaller than Si and S2 when the phases of the Fourier 
components (s r , s^, s z ) of the source density s do not bear 



any relationship with one another. Unlike Si and S2 that 
depend on the absolute values of these (generally complex) 
quantities, S3 sensitively depends on their relative phases. 
In addition, the first two terms are positive for all n, while 
the third term oscillates between positive and negative val- 
ues. Since the measurement of this term by any instrument 
would entail an integration with respect to frequency over a 
non-zero bandwidth, these oscillations will likely make the 
detected value of the third term negligibly smaller than those 
of Si and 5*2- 



3.5 The Rayleigh distance 

The coefficient K r changes from being independent of n to 
decaying as n -1 / 2 when n increases past 7rRp/(r > — f<) 2 
(Table [2). On rearranging this condition for n in terms of 
Rp, we obtain 



Rp 



a 

X 



(6) 



where A = c/v is the wavelength of the electromagnetic ra- 
diation and a = r> — r< is the radial extent of the emitting 
current. In geometrical optics, the parameter c£r is known 
as the Rayleigh distance for an aperture or source of lin- 
ear dimension a; it corresponds roughly to the transition 
between the conditions for Fresnel and Fraunhofer diffrac- 
tion ( |Brooker||2003| ). 

Whilst the Rayleigh distance is most frequently encoun- 
tered in optical or radio experiments when using monochro- 
matic radiation and varying the distance between source and 
observer ( |Brooker]|2003| ), here we cross it by keeping Rp 
fixed and varying the frequency v — nuj/2ix. In the follow- 
ing fits to data, the transition across the Rayleigh distance 
accounts for the steepening of the observed spectra by the 
factor n -1 in the X-ray band and enables us to derive the 
radial extent r> — f < of the emitting plasma that contributes 
towards the subpulse detected at the observation point. 

Having described the general behaviour of equation Q , 
we now embark on the fits of this equation to observational 
data. 



3.6 The Crab pulsar, PSR B0531+21 

We have taken the phase- averaged spectra of the pulsed 
emission that is received from the Crab pulsar in various fre- 



quency bands (the data points shown in Fig. |8| from Lyne 
Graham- Smith (2006). The curve in Fig. [8] represents the 
radiation flux given by equation Q for the following values 
of the parameters: 



1.9 x 10 4 with 



2tt 



30.3 Hz, 



-11/6 
-1/12 



3.2 x 10 b <n< 10 1 



Si oc { n~' /rz 10 1 

v n- 7 / 12 7.9 x 10 16 

$2 S3 

01 01 



< n < 7.9 x 10 1 

< n < 3.2 x 10 2 



(7) 



(8) 



(9) 



for all n, and m — 10 . 

Give n that the Crab puls ar is at the distance Rp c± 6.2 x 
10 21 cm (Taylor et al. 1993) and has a light cylinder with 
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Figure 9. The points show observational data (where available) 
from the spectrum of the Vela pulsar. The curves show the spec- 
tral distribution log(dP n /df2p), predicted by equations j4|, 
and ( |12| ), versus logn and logv for v = noo /(2tt) ~ 11. 2n Hz 
and Q/oo ~ 2.8 x 10 5 . In the model, the recovery of intensity in 
the X-ray band is caused by resonant enhancement due to the 
azimuthal modulation frequency moo / '(2ir) ~ 1.5 x 10 16 Hz. The 
steepening of the gradient of the spectrum by —1 at 3.6 x 10 21 
Hz reflects the transition through the Rayleigh distance. 



log 10 n 

8 10 12 14 16 18 20 22 24 26 
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log 10 v (Hz) 

Figure 10. The points show observational data (where available) 
from the spectrum of the Geminga pulsar. The curves show the 
spectral distribution log(dP n /df2p), predicted by equations |4j, 
( |14| and ( |15| , versus logn and \o%v for v = noo /(2tt) ~ 4.22n Hz 
and fl/oo ~ 5.89 x 10 5 . In the model, the recovery of intensity over 
the optical and X-ray bands is caused by resonant enhancement 
due to the azimuthal modulation frequency moo / (2tt) c± 6.7 x 10 14 
Hz. The steepening of the gradient of the spectrum by —1 at 
4.2 x 10 21 Hz reflects the crossing of the Rayleigh distance. 



the radius c/oo — 1.6 x 10 8 cm, the transition through the 
Rayleigh distance would account for the observed steepening 
of its spectrum at n — 7.9 x 10 16 (Fig. ^ if the radial extent 
r> — f < of the emitting plasma is a fraction 3.9 x 10 -2 of 
the light-cylinder radius. 



3.7 The Vela pulsar, PSR B0833-45 



We have taken the data poi nts shown in Fig.|9]from Thomp-| 
son et al.| fl!999| ) [see also |Lyne fc Graham-Smith fl2006| )]. 
The curve in Fig. [9] represents the radiation flux given by 
equation H for the following values of the parameters: 



- ~ 2.82 x 10 5 

oo 



with — = 11.2 Hz, 
2tt 



(10) 



3.8 The Geminga pulsar, PSR J0633+1746 

We have taken the data points shown in Fig. [10] from 



Thompson et al. (1999) and 



also |Lyne fc Graham- Smith (2006)]. The curve in Fig. 
represents the radiation flux given by equation (EI 



Kargaltsev et al.] ( |2005| ) [see 

To 

for the 



following values of the parameters: 



5.89 x 10° with 



Si oc < 



00 

2^ 



4.22 Hz, 





10 7 


< n < 


1.58 x 10 9 


n 1 / 6 


1.6 x 10 9 


<n< 


1.6 x 10 14 




1.6 x 10 14 


<n< 


3 x 10 16 


n 1 / 4 


3 x 10 16 


< n < 


10 21 




10 21 


< n < 


1.8 x 10 23 


n' 1 


1.8 x 10 23 


< n < 


10 26 , 



(13) 



(14) 



n 



Si oc < 



-1/2 
3/2 



n 

n 
n 



S2 

Si 



io 7 , 



•3/10 
•1 

Ss 



2.5 x 10 7 
1.3 x 10 9 

1.3 x 10 15 
4 x 10 16 

3.2 x 10 20 
10 23 

< 1 



< n < 

< n < 
<n< 
<n< 

< n < 
<n< 



1.3 x 10 9 
1.3 x 10 15 
4 x 10 16 
3.2 x 10 20 
10 23 

1.6 x 10 26 , 



(11) 



(12) 



for all n, and m = 1.3 x 10 15 . 

Give n that the Vela puls ar is at the distance Rp ~ 1.5 x 
10 21 cm (Tay lor et al.||l993 ) and has a light cylinder with 
the radius c/oo — 4.3 x 10 cm, the transition through the 
Rayleigh distance would account for the observed steepening 
of its spectrum at n — 3.2 x 10 20 (Fig. [9]) if the radial extent 
r> — f< of the emitting plasma is a fraction 1.9 x 10 -4 of 
the light-cylinder radius. 



& 
Si 





10 7 



1.6 x 10 14 < n < 3 x 10 16 
otherwise, 



(15) 



S3 < Si for all n, and m = 1.6 X 10 14 . 

Given that the Geminga pu lsar is at the distance Rp ~ 
4.6 x 10 20 cm ( |Taylor et aLl|l993| ) and has a light cylinder 
with the radius c/oo = 1.1 x 10 cm, crossing the Rayleigh 
distance would account fo r th e observed steepening of its 
spectrum at n — 10 21 (Fig. 10) if the radial extent f> — f< 



of the emitting plasma is a fraction 3.6 x 10 5 of the light- 
cylinder radius. 



3.9 PSR B0656+14 

We have taken the data points shown in Fig. 



Shibanov et al. (2006), De Luca et al 



galtsev et al. (2007). The curve in Fig 



2005), and Kara- 



m 

id[I 



from 



11 represents the 
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Figure 11. The points show observational data (where available) 
from the spectrum of B0656+14. The curves show the spectral 
distribution log(dP n /df2p), predicted by equations j4|, ( |17| and 
( |18[ ), versus logn and logz/ for v = nuj /(2tt) ~ 2.6n Hz and 
Q/uj ~ 5.89 x 10 5 . In the model, the recovery of intensity over the 
optical and X-ray bands is caused by resonant enhancement due 
to the azimuthal modulation frequency moo /(2tt) ~ 5.2 x 10 14 
Hz. The steepening of the gradient of the spectrum by —1 at 
1.8 x 10 19 Hz corresponds to crossing the Rayleigh distance. 
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Figure 12. The points show observational data (where avail- 
able) from the spectrum of B1055-52. The cuves show the spec- 
tral distribution log(dP n /df2p), predicted by equations j4|, ( |20| 
and \2\\ . versus logn and log v for v = nuo / (2tt) ~ 5.07n Hz 
and Q/uj ~ 6.31 x 10 5 . In the model, the recovery of intensity in 
the optical band is caused by resonant enhancement due to the 
azimuthal modulation frequency moo / (2ir) ~ 6.6 x 10 14 Hz. The 
steepening of the gradient of the spectrum by —1 at 4.5 x 10 20 Hz 
is due to the crossing of the Rayleigh distance. 



radiation flux given by equation Q for the following values 
of the parameters: 



- ~ 7.94 x 10 5 with 



2tt 



2.6 Hz, 



(16) 



Si oc 



^3/8 
1/10 



-1/2 



6.3 x 10 7 
1.3 x 10 14 
3.7 x 10 16 
8.9 x 10 19 



<n< 1.3 x 10 14 

<n< 3.7 x 10 16 

<n< 8.9 x 10 

<n< 10 24 



19 



(20) 



2~ 1/2 10 7 < n < 3.2 x 10 12 

i~ 3/8 3.2 x 10 12 < n < 2 x 10 14 

Si oc { n 3/8 2 x 10 14 < n < 7.6 x 10 16 

i° 7.6 x 10 16 < n < 7.1 x 10 18 

i- 1/2 7.1 x 10 18 < n < 1.58 x 10 23 , 



& _ f 10 7 < n< 7.6 x 10 16 

Si~ \ 1.78 x 10 6 otherwise, 



(17) 



(18) 



Ss < Si for all n, and m = 2 x 10 14 . 

Give n that PSR B0656+ 14 is at the distance R P ~ 2.3 x 
10 21 cm (Tay lor et gtl.||l993] ) and has a light cylinder with 
the radius c/uj — 1.8 x 10 cm, the transition through the 
Rayleigh distance would account for the observed steepening 



of its spectrum at n — 7.1 x 10 (Fig. 11 ) if the radial extent 



r> — f< of the emitting plasma is a fraction 7.5 x 10 
the light-cylinder radius. 



of 



3.10 PSR B1055-52 

We have taken the data points shown in Fig. [12 from 
Thompson et al.|(|1999|) and |McCulloch fe Hamilton| ( 1976]) 



see also Lyne & Graham-Smith (2006 )]. The curve in Fig. 12 
represents the radiation flux given by equation Q for the 
following values of the parameters: 



6.31 x 10 & with 



2tt 



5.07 Hz, 



(19) 



S2 _ \ 1.3 x 10 14 < n< 3.7 x 10 1 

Si I 4.22 x 10 5 otherwise, 



(21) 



: 4.7 x 



S 3 < Si for all ra, and m = 1.3 x 10 14 . 

Give n that PSR B 1055-5 2 is at the distance Rp 
10 21 cm ( [Taylor et al. 1993) and has a light cylinder with 
the radius c/cu = 9.4 x 10 cm, the transition through the 
Rayleigh distance would account for the observed steepening 



of its spectrum at n = 8.9 x 10 (Fig. 12 ) if the radial extent 



r> — r< of the emitting plasma is a fraction 4.2 x 10 
the light- cylinder radius. 



of 



3.11 PSR B1509-58 

We have taken the data points shown in Fig. |13| from 
Thomps on et ah] ( |1999| [see also Lyne & Graham-Smith 



(2006)]. The curve in Fig. 13 represents the radiation flux 



given by equation Q for the following values of the param- 
eters: 



5.62 x 10 & with 



2^ 



6.6 Hz, 




5 x 10 7 <n< 2.5 x 10 15 
2.5 x 10 15 < n < 5 x 10 18 
5xl0 18 <n< 1.6 xlO 26 , 



<T<1, f «1 



(22) 

(23) 
(24) 
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Figure 13. The points show observational data (where available) 
from the spectrum of B 1509-58. The curves show the spectral 
distribution log(dP n /df2p), predicted by equations j4|, ( |23| and 
( |24| , versus logn and logz/ for v = noo /(27r) ~ 6.6n Hz and 
fi/u; ~ 5.62 x 10 5 . The steepening of the gradient of the spectrum 
by —1 at 3.3 x 10 19 Hz is caused by crossing the Rayleigh distance. 



for all n. (The existing data is insufficient to pinpoint the 
value of m for this pulsar.) 

Give n that PSR B1509-5 8 is at the distance R P ~ 1.3 x 
10 22 cm ( Taylor et al.||l993 ) and has a light cylinder with 
the radius c/oo = 7.2 x 10 s cm, the transition through the 
Rayleigh distance would account for the observed steepening 
of its spectrum at n = 5 x 10 18 (Fig. 13) if the radial extent 
r> — f < of the emitting plasma is a fraction 3.4 x 10 -3 of 
the light-cylinder radius. 



3.12 PSR B1706-44 

We have taken the data points shown in Fig. |14| from 
|Thompson et ah] ( |1999| [see also |Lyne &; Graham-Smith] 
(2006)]. The curve in Fig. 14 represents the radiation flux 



given by equation Q for the following values of the param- 
eters: 



- ~ 3.71 x 10 5 

00 



with 



2tt 



Si oc < 



-1/2 



,1/2 



3.2 x 10 7 
3.2 x 10 14 



9.76 Hz, 



< n < 3.2 x 10 14 

< n < 8.3 x 10 15 



-1/10 



^22 



(25) 



(26) 



Figure 14. The points show observational data (where available) 
from the spectrum of B 1706-44. The curves show the spectral 
distribution log(dP n /df2p), predicted by equations j4|, \26\ and 
( |27| , versus logn and \ogu for v = noo / (2tt) ~ 9.76n Hz and 
fl/oo ~ 3.71 x 10 5 . In the model, the recovery of intensity in 
the optical band is caused by resonant enhancement due to the 
azimuthal modulation frequency moo / (2ir) ~ 3.1 x 10 15 Hz. The 
transit through the Rayleigh distance accounts for the steepening 
of the gradient of the spectrum by —1 at 1.9 x 10 23 Hz. 



3.13 PSR B1929+10 



We have taken the data points shown in Fig.[l5]from Becker 
et al. (2006). The curve in Fig. 15 represents the radiation 
flux given by equation Q for the following values of the 
parameters: 



- - 6.3 x 10 5 with — = 4.41 Hz, 



(28) 



51 oc n" 1/2 , = 3.16xl0 6 , 2xl0 7 < n< 1.6 x 10 14 , (29) 

oi 

5 2 oc n 3/8 , I 1 < 1, 1.6 x 10 14 < n< 2.8 x 10 17 (30) 

02 

53 is negligible for all n, and m = 1.6 x 10 14 . 

Given tha t PSR B1929+10 i s at the distance R P ~ 
5.2 x 10 20 cm ( |Taylor et a]T||l993| | and has a light cylinder 
with the radius c/00 — 1.1 x 10 cm, the absence of a sudden 



< n < 



10 



27 



steepening in the gradient of its spectrum (Fig. 15) means 
that the transition through the Rayleigh distance must occur 
at harmonic numbers larger than n — 2.8 x 10 17 , and so the 
radial extent f> — f< of the emitting plasma must be smaller 
than a fraction 2.3 x 10 -3 of the light-cylinder radius. 



S 2 
01 



01 



(27) 



for all n, and m = 3.2 x 10 14 . 

Give n that PSR B 1706-4 4 is at the distance R P ~ 5.6 x 
10 21 cm (Tay lor et al.||l993 ) and has a light cylinder with 
the radius c/00 — 4.9 x 10 cm, the transition through the 
Rayleigh distance would account for the observed steepening 



of its spectrum at n = 2 x 10 (Fig. 14) if the radial extent 



r> — f< of the emitting plasma is a fraction 4.2 x 10 
the light-cylinder radius. 



of 



3.14 PSR B1951+32 



We have taken the data points shown in Fig. [16] from 
|Thompson et al.] ( |1999| [see also |Lyne &; Graham-Smith] 



(2006)]. The curve in Fig. 16 represents the radiation flux 
given by equation Q for the following values of the param- 
eters: 



- ~ 1.41 x 10 5 

00 



with 



2tt 



25.3 Hz, 



(31) 
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Figure 15. The points show observational data (where avail- 
able) from the spectrum of B1929+10. The curves show spectral 
distribution log(dP n /df2p), predicted by equations j4|, ( |29| and 
( |30| ), versus logn and logz/ for v = nuj /(2tt) ~ 4.41n Hz and 
~ 6.3 x 10 5 . In the model, the recovery of intensity in the 
optical band is caused by resonant enhancement due to the az- 
imuthal modulation frequency muo /(27r) ~ 7 x 10 14 Hz. 



Table 3. Summary of Q values used in fitting pulsar spectra. The 
parameter / p is the derived plasma frequency, where we assume 
that f p = Q/2tt. N e is the electron density in the emitting region 
of the pulsar's atmosphere calculated from the plasma frequency 
f p , using (27t/ p ) 2 = 47rA/ r e e 2 /m e , where e is the electron's charge 
and m e is its mass. 



Pulsar 


n 

UJ 


UJ 

2tt 

(Hz) 


fp 
(MHz) 


(10 4 cm- 3 ) 


PSR B0531+21 


19000 


30.3 


0.58 


0.4 


PSR B0833-45 
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11.2 


3.2 


12 


PSR J0633+1746 


589000 


4.22 


2.5 


7.7 


PSR B0656+14 


794000 


2.6 


2.1 


5.3 


PSR B1055-52 


631000 


5.07 


3.2 


13 


PSR B1509-58 


562000 


6.6 


3.7 


17 


PSR B1706-44 


371000 


9.76 


3.6 


16 


PSR B1929+10 


630000 


4.41 


2.8 


9.6 


PSR B1951+32 


141000 
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3.6 


16 
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Figure 16. The points show observational data (where available) 
from the spectrum of B1951+32. The curves show the spectral 
distribution log(dP n /df2p), predicted by equations j4|, ( |32| and 
( |33| , versus logn and logz/ for v = nuj /(2ir) ~ 25. 3n Hz and 
Q/uj ~ 1.41 x 10 5 . In the model, the recovery of intensity in 
the optical band is caused by resonant enhancement due to the 
azimuthal modulation frequency muj/(27r) ~ 1.4 x 10 14 Hz. The 
crossing of the Rayleigh distance accounts for the steepening of 
the gradient of the spectrum by —1 at 10 20 Hz. 



( n- 1 / 2 


10 7 


< n < 


5.6 x 10 12 




5.6 x 10 12 


<n< 


5.6 x 10 15 


< n 5 / 12 


5.6 x 10 15 


< n < 


4 x 10 18 




4 x 10 18 


< n < 


10 23 


k n- 2 / 3 


10 23 


<n< 


5 x 10 24 , 



< 1, 



s 3 



< 1 



5*2 S2 

for all n, and m = 5.6 x 10 12 . 



(32) 



(33) 



Given tha t PSR B01951+32 is at the distance R P ~ 
7.7 x 10 21 cm ( |Taylor et aT]|l993| ) and has a light cylinder 
with the radius c/uj = 1.9 x 10 cm, the transit through the 
Rayleigh distance would account for the observed steepening 
of its spectrum at n = 4 x 10 18 (Fig 
r> — r< of the emitting plasma is < 
the light- cylinder radius. 



16 ) if the radial extent 
fraction 5.6 x 10~ 3 of 



4 DISCUSSION 

4.1 Fitted values of Q and inferred plasma 
frequencies 

The fitted values of Q/cu taken from Section [3] are sum- 
marized in Table [3] where we have also converted the 
Q/2tt values into Hz using the known rotation frequen- 
cies of the pulsars. It will be remembered that Q/2tt 
represents a characteristic frequency that modulates the 
source [see equation (|3|] in the emitting region; assuming 
that the pulsar magnetosphere consists of a dilute neutral 
plasma ( |Spitkovsky||2006| ) , it is reasonable to make the at- 
tribution Q/2-7T = / p , where f p is the (electron) plasma fre- 
quency in the emitting part of the magnetosphere. The fifth 
column of Table|3]gives inferred electron densities iV e derived 
using (27t/ p ) 2 = 47riV e e 2 /rae, where e is the electron's charge 
and m e is its mass ( |Lyne &; Graham-Smith 2006). Note 
that virtually all of the electron densities are of the order of 
10 5 cm~ 3 , with only the Crab's value (iV e = 4x 10 3 cm~ 3 ) be- 
ing somewhat lower [^Nevertheless, all of these iV e are com- 
pletely consistent with the densities expected in a conven- 
tional pulsar magnetosphere ( |Lyne &; Graham-Smith|2006| ). 
Another way of visualizing the rather similar values of 



1 The lower electron density around the Crab could result from 
its very high optical luminosity ( |Lyne fe Cr aham-Smith 2006), 
which might act to "blow" electrons out of the magnetosphere. 
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Q found in Sectio n |3| for all of the pulsars is to plot 0,/cu 
versus cu/2n (Fig. 17); the values of Q/uj all lie close to the 
curve n/cu = 2.56 x 10 6 (2tt/u;). 

Given that positions of the last peaks of the Airy func- 
tions appearing in equation Q and hence the salient fea- 
tures of the pulsar spectrum scale as (Q/cj) 3 , the rather 
uniform values of imply that slower pulsars (those with 
smaller values of co) will have observed spectral intensities 
weighted towards higher frequencies. By contrast, millisec- 
ond pulsars (with large u) should have emission concen- 
trated at lower frequencies than the above pulsars. These 
predictions seem to be borne out both by the radio-quiet 
pulsar in the supernova remnant CTA1 (uj/2tt ~ 3.1 Hz), 
which has emission peaked in the gamma-ray end of the 
spectrum ( Abdo et al. 2008), and by millisecond pulsars such 
as J1748-22446ad (oj/2tt ~ 716 Hz), B1937+21 (uj/2tv ~ 642 
Hz) and B 1957+20 (uj /2tt ~ 622 Hz) which show no emission 
in the high-frequency range but generally strong and bright 
pulses at radio frequencies ( Soglasnov et al.[20 04; Hessels et 



aLl[2006l |Huang fe Becker||2007 ) 



4.2 The frequency muj /2tv and inferred magnetic 
fields 

Table[4]shows the values of m derived from the frequencies at 
which the observed spectra recover their intensity as a result 
of resonance enhancement (see Section [3]) . The frequency 
mcj/27T represents the frequency of azimuthal fluctuations of 
the polarization-current distribution pattern. These spatial 
fluctuation frequencies (~ 10 13 — 10 15 Hz) could well arise 
from the cyclotron resonance of electrons. The cyclotron (an- 
gular) frequency co c in a field H is defined as oo c = eH/cm e , 
where m e is the mass of the electron and e is its charge (cgs 
units) ( |Lyne &; Graham-Smith|2 006); hence, we can use the 
mu values to infer magnetic fields of around 10 7 — 10 9 G 
in the emitting region (Table |4j). Such fields are completely 
consistent with those expected in a conventional pulsar mag- 



netosphere (Lyne & Graham-Smith 2006). 



4.3 The size of the emitting region 

The set of frequencies at which the slopes of the observed 
spectra suddently steepen by —1 mostly lie in the range 
10 18 — 10 21 Hz. This implies that the radial extent of the 
emitting polarization current is a fraction 10 -2 — 10 -5 of the 
light-cylinder radius (Section^. Recalling that the subbeam 
that arises from the emitting part of the source in the radial 
interval f< ^ r ^ f> (where r = ru/cj is detectable within 
the polar interval arccos(l/f<) ^ \6p — tt/2\ ^ arccos(l/r>) 
(see Section 2.3.1), it follows that the polar widths A0p ~ 
sin Op tan 0p{r > — f<) of the subbeams that are detected 
outside the plane of rotation (Op ^ tt/2) are of the order of 
10 -2 — 10 ~ 5 radians. [Note that the overall radiation beam 
consists of a superposition of such narrow subbeams (Sec- 
tion EO.] 



4.4 



Frequency dependence of Fourier components 
of the source density s 



Table 4. Summary of m values used in fitting pulsar spectra. 
The final column tabulates the magnetic fields deduced on the 
assumption that mw /2tt represents the cyclotron resonance of 
electrons in the emitting region of the magnet osphere. These were 
derived using the standard equation for the cyclotron (angular) 
frequency uj c : muj = uj c = eH/cm e , where e is the electron's 
charge and m e is its mass. 



Pulsar 




m 




muj 
2tt 


H 












(Hz) 


(THz) 


(G) 






PSR B0531+21 


10 12 


30.3 


30 


1.1 


x 


10 7 


PSR B0833-45 


1.3 


x 10 15 


11.2 


15000 


5.2 


X 


10 9 


PSR J0633+1746 


1.6 


x 10 14 


4.22 


680 


2.4 


X 


10 8 


PSR B0656+14 


2 x 10 14 


2.6 


520 


1.9 


X 


10 8 


PSR B1055-52 


1.3 


x 10 14 


5.07 


660 


2.4 


X 


10 8 


PSR B1509-58 






6.6 










PSR B1706-44 


3.2 


x 10 14 


9.76 


3100 


1.1 


X 


10 9 


PSR B1929+10 


1.6 


x 10 14 


4.41 


710 


2.5 


X 


10 8 


PSR B1951+32 


5.6 


x 10 12 


25.3 


140 


5.1 


X 


10 7 




The dependence of the Fourier components s r ,^, z of the 
source densities s r ,ip, z on frequency (or on harmonic number 



15 20 

uj/2tt (Hz) 



Figure 17. The points designate the values of the two parameters 
Q/uj and uj/2tt for the 9 pulsars considered in Section [3] and the 
curve Q/u; = 2.56 x 10 6 (27r/o;) represents the best fit to these 
points. 



n) is determined by the fluctuations (within the emitting 
region) of s r ,^>, z with the coordinate z [see equation (A27)]. 
The small-scale spatial variations of s in the direction paral- 
lel to the rotation axis that are implied by the observational 
data (see Section [3| may be identified with the sharp gra- 
dients in the distribution of the source density at a current 
sheet. The two-dimensional current sheets predicted by nu- 
merical models of MHD pulsar magnetospheres extend be- 
yond the light cylinder. Not only do they have distribution 
patterns that rotate superluminally (Spitkovsky 2006), but 
they also vary discontinuously, and so have fluctuations that 
occur over wide ranges of length scales and amplitudes. 
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5 SUMMARY AND CONCLUSION 

We have compared the observed spectra of 9 pulsars for 
which multiwavelength data are available with the spectrum 
of the radiation generated by an extended source in superlu- 
minal rotation and shown that the entire spectrum of each 
pulsar can be accounted for quantitatively in terms of this 
single emission process. The results reported here are model- 
independent in that the only global property of the magne- 
tospheric structure invoked is its quasi-steady time depen- 
dence: the cylindrical components j r ,ip,z 0", z; t) of the den- 
sity of the magnetospheric polarization current depend on 
ip only in the combination tp — uot. This property follows un- 
ambiguously from the observational data and implies that a 
current distribution with a superluminally rotating pattern 
at a radius r > c/co is responsible for the unique features 
of pulsar emission such as the pulses' extreme brightness 
temperature, temporal width, source dimension and peak 
spectral frequency as well as the average pulses' polariza- 
tion properties (their occurrence as concurrent 'orthogonal' 
modes with swinging position angles and with nearly 100 
per cent linear or circular polarization). 

The curves in Figs. |&ffT6| demonstrate that the universal 
features of pulsar spectra (their breadth over 15-18 orders of 
magnitude of frequency, the recovery of their amplitudes in 
the optical band, their bell-shaped peaks in ultraviolet, and 
the sudden steepening of their gradient by —1 in the X-ray 
band) are all described by the properties of the Green func- 
tion for the problem, and so are consequences solely of the 
accelerated superluminal motion of the distribution pattern 
of their source. Only the finer details of these spectra depend 
on the structure of the pulsar magnet osphere, or more specif- 
ically, on the structure of the current sheets just outside the 



light cylinder of the pulsar magnetosphere (ppitkovsky 2006 ) 



(i.e. on the variable 'weather' in the pulsar magnetosphere, 
as opposed to its stable 'climate'). 

The increase in the spacing between successive peaks of 
the oscillating spectra shown in Figs. |8lfl6l with frequency, 
which has already been observed in the radio emission from 
the Crab pulsar ( |Hankins &; Eilek||2007 ), also only depends 
on the properties of the Airy functions in equation Q; it 
is another universal feature of the pulsar emission predicted 
by the present model that can be tested, in principle, over 
any frequency band [see Fig. [7] and Ardavan et al. (2008a)]. 
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APPENDIX A: MATHEMATICAL 
DESCRIPTION OF THE RADIATED FIELD 

A general mathematical treatment of spectrum of the radi- 
ation field that is generated by the source distribution de- 
scribed in equation ([3} has already been given in Ardavan et 
al. ( 2003 ). Our purpose in this appendix is to make the phys- 



ical content of the previous analysis more transparent by 
presenting its results from an alternative point of view. We 
begin by calculating the spectrum of the radiation emitted 
by one of the constituent rings of source by first deriv- 
ing the spectrum of the Lienard-Wiechert potential due to a 
superluminally rotating point source, and then superposing 
the potentials of the point-like line elements that make up 
the ring. Once we have clarified the origin of the higher fre- 
quencies in the emission from this simpler ring current, we 
will further superpose the fields of the rings that make up 
the volume-distributed source ^ to arrive at the spectrum 
on which the modelling of the observed spectra considered 
in this paper is based. 



Al Radiation field of a constituent ring of the 
source 

The trajectory of a uniformly rotating point source, e.g. the 
source of synchrotron radiation, is described by 



<p = + ujt, 



(Al) 



where the subscript denotes the constant values of the 
cylindrical polar coordinates (r, ip, z) of the point source at 
time t = 0, and co is constant. (Note that, to describe re- 
peated rotations, ip in this equation increases indefinitely as 
the time t does.) If this source has a strength s that oscil- 
lates with the frequency Q in time, then the density of the 
polarization current it carries can be written as 



j = six} exp(iQ£)(*)(r — ro)6(z — zo)S(cp — (fo — ufje^, (A2) 

where 6 is the Dirac delta function, and e^ is the base vector 
associated with the azimuthal coordinate ip. 



The vector potential due to this oscillating, rotating 
point source is given by 



point 



-U 



rdrdcpdzdt j 



S(t P - t-R/c) 
R 



/ d£exp(i£2£) — ^^e^o, (A3) 

c J 



in which 
Ro = R\ 
with 



Ro 



^ip \r=r ,(p=(po+ujt,z=zo 5 



(A4) 



R= [(z P - zf +r P 2 + r 2 



■ 2rrp cos((pp — ip)] 



1/2 



and 



: s'm((p P - (p)e rp + cos((pp - (p)e^ 



(A5) 



(A6) 



where (e rp , e^p , e Zp ) are the constant base vectors associ- 
ated with the coordinates (rp,(pp, zp) at the observation 
point P, and c is the speed of light in vacuo. Evaluation of 
the remaining integral with respect to t yields 



> point 



srooo 



E 



exp(i^t) 



\Ro - (rorpu/c) sin((p P — ipo — u)t)\ 



x [sm(cpp 
+ cos((p P 



■ out)e rp 



where the retarded times tj are the solutions of 



t P = t + R /c 



(A7) 



(A8) 



[see, e.g. |Jackson (1999)]. Depending on the values of 
(ro, (po, zo), there are either one or three tjS when the speed 
rooo of the source exceeds c moderately [see Fig. [l] and Ar- 



davan|fll998| )]. 

Taking the curl of A pomt , and noting that the depen- 
dence of this potential on the coordinates of a far-field ob- 
server arises primarily from the dependence of the retarded 
times tj on these coordinates, we obtain the following ex- 
pression for the magnetic field of the radiation emitted by 
the above point source: 



>point 



srooo ^ 

'~c^~ 2^ 



dt 



1 - r fp sin((^p - ip — ut)/Ro 
exp(iQt)p 



Ro - r rp sm((p P — (po — ui) 



(A9) 



with 



p = cos#p sm(cpp — cpo — ujt)e\\ +cos((^p — ipo —ut)e±, (A10) 

where ey = e^p (which is parallel to the plane of rotation) 
and e^ = nxey comprise a pair of unit vectors normal to 
the radiation direction 



n = sin#pe rp + cos 6pe 



(All) 



and ro = rou/c, rp — rpu/c, and Ro = Rou/c. This is 
the familiar Lienard-Wiechert field, except that it receives 
contributions from more than one retarded time tj if ro > 1 
and so the source moves faster than light (see Section^. 

Next, let us consider a rotating ring r = ro, z — zo, 
whose strength varies spatially, as cos(m(^o), in addition to 
oscillating temporally (m is an integer). The field due to 
such a ring can be obtained from the field B pomt of the 
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above point source by multiplying B pomt by cos(mupo) and 
superposing the contributions of the elements that make up 
the ring: elements whose azimuthal positions cpo at t — 
cover a full circle only once, and so are labelled uniquely 
by their initial coordinates over the interval — ty < ipo ^ tt. 
That is to say, 



/7I 
-1 



B 



point 



cos(m<po)dv? , 



(A12) 



with a range of integration that is limited to 2tt. [Note that 
the variable opo is the same as the variable up appearing in 
equation 



For ro < 1, the denominators in equation (A9) nowhere 
vanish, and so B pomt is a regular and periodic function of 
(po + ootp: according to equation (|A8|), tj — tp are periodic 



functions of ipo+wtp with the period 2tt/oo. In this sublumi- 
nal case, the expression on the right-hand side of equation 



(A12) constitutes the coefficient of a single term in the co- 
sine series for B pomt . Thus the Fourier expansion of B rmg 
consists of only two sinusiodal functions of opo + cotp with 
the frequencies fi ± moo. For fo > 1, on the other hand, 



the denominators in equation ( A9 ) vanish for a source point 
(ro,(po,zo) that approaches the observer with the speed of 
light and zero acceleration at the retarded time: the ex- 
pression appearing in the first denominator i n this equation 
equals 1 + c~ 1 dRo/dt. Being non-integrable (Ardavan et al. 
2004a), this singularity makes it impossible to Fourier an- 



alyze the representation ( A9 ) of the Lienard- Wiechert field 



B pomt directly to obtain the spectral distribution of this 
field. 

The way to handle the singularity of B pomt in the su- 
perluminal case is to go back one step and instead expand 



the integrand in equation (A3) in a Fourier series with re- 
spect to ipo. (We must use a Fourier series, rather than a 
Fourier integral, because the values of opo over which this 
function is defined lie in an interval of finite length.) The 
result is 



point 



STqOO 

2ttc 



oo ~ 

exp(ino;£p) / dt Ro -1 



x exp{ — i [noo (t + Ro/c) — fit] je^ , 



(A13) 



as can be readily seen by replacing the delta function in 



equation (A3) by its Fourier-series representation. Taking 
the curl of this, we then obtain 



jpoint 



27ric 2 ^ 



nexp(ino;tp) J dt Ro 1 



x exp{— i[nuj(t + Ro/c) — fit]} 
x [cos Op sin((pp — cpo — ujt)e\\ 
+ cos((^p - <po - oot)e±]. 



(A14) 



For a radiation frequency noo/2ix that appreciably exceeds 
the rotation frequency, the asymptotic value of the integral 
in the above expression receives contributions solely from 
the stationary points of the phase, ct + Ro(t), of the rapidly 
oscillating exponential in its integrand. 

This phase is stationary, in the course of each rotation, 
at the following two retarded times at which the source ap- 
proaches the observer along the radiation direction with the 
speed c: 



ujt± = UDp 



3tt , 

(fo + — ± arccos 



1 )+2/ctt, (A15) 
\r sin Up J 



K ro sine 

where (Rp,0p, opp) (with Rp ^> c/oo) denote the spheri- 
cal polar coordinates of the observation point and k is an 
integer. These stationary points coalesce if the observer is lo- 
cated at Op — arcsin(l/fo), i.e. when the source approaches 
the observer not only with the speed c but also with zero 
acceleration [see Fig. [I] and Ardavan et al. ( 2007)]. In this 
case, the rapid oscillations of the exponential in equation 



(A14) result in the destructive interference of contributions 



from all emission times, except those made during a short 
interval centred at t = t c , where uot c — Lpp — (/?o + 37r/2 + 2/c7r. 

For n ^> 1 and Op — arcsin(l/fo), we can therefore 
obtain the leading term in the asymptotic expansion of the 
integral in equation (A14) by applying the principle of sta- 



tionary phase: by approximating both the phase ct + Ro(t) 
and the amplitude of the rapidly oscillating exponential with 
the dominant terms in their Taylor expansions about t — t c 



and replacing the limits of integration with zboo (Borovikov 



1994). Evaluation of the resulting integral with the aid of 
equation (10.4.32) of |Abramowitz Stegun| ( |1970| ) thus 
yields 



jpoint 



ic 2 R P 



E 

|n|>l 



2/3 



exp{i[n(o;tp — Rp + zo cosOp) 
+ (n - Q/u)((po -tpp- 3tt/2)]} 
' 2 \ 1/3 fi" 



{ 



cos #pAi 



(!) 



+i I 



where Ai and Ai ; stand for the Airy function and the deriva- 
tive of the Airy function with respect to its argument, re- 
spectively, and Rp = Rpoo/c and zq = zquj / c. 



Amplitudes of the individual terms in equation (A16) 
have the dependence n -1 / 4 on n for n 1, so that this 
series yields a divergent value for B pomt on the cusp locus 
Op — arcsin(l/ro), as did equation (A14). However, these 
individual Fourier components are not periodic functions of 
(po when fi differs from an integral multiple of a;. As a re- 



sult, when we insert equation (A16) in equation (A12) and 



carry out the integration with respect to <^o, we obtain an 
expression that contains all values of (n ^> 1): 



Ed) 



2/3 



ic 2 R P 

|n|»l 

exp{i[n(u;tp — Rp + zq cosOp) 
-(n-fi/cj)((^P+37r/2)]} 

'2y/ 3 fi~ 



x < cos Op Ai 



where 



sin[(n + m — Q/lj)tt] 

n + rn — fi/o; 
^sin[(n — m — fi/o;)7r] 
n — m — fi/o; 



(A18) 
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is the factor resulting from the integration over ipo. This 
series has an infinite number of non-zero terms because, on 
the one hand, the variable ipo that labels the elements of 
the ring can only range over an interval of length 27T, and 
on the other hand, the function B pomt that we are Fourier 
analyzing has different values at opposite ends of such an 



interval [see Section 2 of Ardavan et al. (2003)]. 

When Q/u; is an integer, only the two terms n — =bm + 
Q/uj of the series in equation (A17) are non-zero: the nu- 
merators and denominators of the two fractions in equation 



(A18) simultaneously vanish approaching the limits tt. But 



when Q/co is nonintegral, this series diverges, as expected 
from the fact that the singul arity of the time-domain ex- 
pression for B pomt in equation (A9) is non-integrable. There 
is a radical difference between the series in equation (A17) 



and its subluminal counterpart. Not only is it not possible to 
derive the frequency-domain expression for B rmg by directly 
Fourier analyzing the time-domain expression for this field 
[which is divergent when the distance between the observer 
and the source decreases with the speed dRo/dt = — c], 
but, correspondingly, the transformation t — > t + Ro(t)/c 
that takes one from the frequency- domain expression to the 
time-domain expression has a vanishing Jacobian and so is 
mathematically impermissible in the superluminal case. a25 
On the one hand, each source element ipo of the ring 
makes its contribution towards the observed field at the dis- 
crete set of retarded times cut = ipp — ipo + 3tv/2 + 2kit at 
which it approaches the observer with the speed of light and 
zero acceleration, i.e. periodically with the period 2tv/co. On 
the other hand, temporal modulations of the density of the 
source occur on a time scale 2tt/Q that is incommensurable 
with the period 2n/cv. Because the emission time of the ele- 
ment labelled by ipo is fixed (via cot = ipp — ipo + 3tv/2 + 2kn) 
by its initial azimuthal position (ipo), the temporal [exp(iQt)] 
and spatial [cos(mipo)] modulations of this source effectively 
combine into a single variation [oc exp(iQipo/uj) cos(mipo)]. 
It is the Fourier decomposition of this variation with ipo, 
or equivalently with the retarded time t, that together with 
the incommensurablity of the two periods (2tt/lv and 2ir/Q), 
and the finiteness of the domain of definition of cpo , results in 
a spectrum containing all frequencies [see Fig. 5 of | Ardavan] 
et al.|(|2003|]. 



A2 Radiation field of the entire volume of the 
source 

In this Subsection, we use a similar procedure to find the 
field generated by the volume-distributed source described 
in equation ^ for which 



j = dP/dt = |is(r, z) J^±(^ ± mu) exp[i(Qt =p mipo) 



+{m 



,0 



(A19) 



where the symbol {m — > —m, Q — > —£7} designates a term 
exactly like the one preceding it but in which m and Q are 
everywhere replaced by — m and — Q, respectively. Here, we 
have replaced ip by ipo to bring out the connection with the 
ring source considered in Subsection |A1[ 



Taking the curl of the first member of equation (A3) 
and discarding terms of order R~ 2 , we obtain the following 



generalization of equation (A14): 



— / rdrdipdzdt — ^^nxj 



R 



2tvc 2 



oo 

nexp(ino;tp) / rdrdifodzdtRo -1 



x exp[— in{ujt + i2o)]nxj, 



(A20) 



where 



Ro = [(zp-z) +r P +r -2rr P cos((p P -(p -ut)] 1 , (A21) 

as before. Note that the change of integration variable from 
(p to (fo = (p — Lot in the second member of this equation, that 
is made to emphasize the correspondence between equations 



( A14) and ( A20 ), has a non-zero Jacobian and so is perfectly 
permissible. 

For n ^> 1 and a far-field observer outside the plane of 
rotation, i.e. Op / 7r/2, the phase t + Ro(t) / c of the rapidly 
oscillating exponential in equation ( A20 ) is stationary at t = 
t Cl r = esc Op, simultaneously, where cvt c = tpp — (^o + 37r/2 + 
2/c7r as in Subsection |A1| Performing the integration with 
respect to these two variables by the method of stationary 
phase, we obtain the following expression for the electric 
field E = — hxB of the radiation: 



E = ^ I E + 2 ^ En exp(-m^p) 

I n=l 

in which 

m — Q/uj / n \ 2 / 3 



(A22) 
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with 

K r ~ 
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f> — f<, n < 7rRp/(f > - r<) 



(27iR P /n) 1/2 exp(-i7r/4), n > ttR p / '(f>-f<) 2 , (A25) 
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(A26) 



(A27) 



Here, stands for the real part of the quantity inside brack- 
ets, s r ,ip, z are the cylindrical components of s, and f< and f> 
denote the radial boundaries of the superluminal part of the 



source that contributes towards the radiation at P ( Ardavan 
eTaLl[2003 ) 



As in the case of the ring source, the fact that the ra- 
diation observed at the time tp is determined almost exclu- 
sively by the state of the emitting current at the retarded 
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time t — t c turns the temporal modulation exp(iQt) of 
the source into a spatial modulation exp(— [Q(po/uu). Conse- 
quently, the remaining integral in equation ( A23 ) again has 
a non-zero value for all n when Q/uj differs from an integer. 
Insertion of the resulting expression for E n in dP n /dQp = 
cR P 2 : |E n | 2 /(27r) now yields the radiated power per harmonic 
per unit solid angle: 



dP n 

where 



27TUJ 2 



esc 2 OplKr^K^ 2 
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+ — 

li+ n 



(A28) 



(A29) 



and ji± — zbm + Q/uj. Since the Airy functions in equation 
(A26) decay exponentially with n when their argument is 
positive ( Abramowitz &; Stegun|1970 ), the terms associated 
with — Q < are ignored in the above expression. |^] 

Evaluating the absolute value of the complex vector V 
in equation ( A26 ) , we obtain 
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(A30) 



in which 
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^| + \s r \" COS 2 0p) 1/2 , 

^'i^i^o [s* cos Op (sp cos Op 
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(A32) 



(A33) 



Here, and the superscript star denote the imaginary part 
and the conjugate of a complex variable, respectively. 



2 This expression for dP n /dflp slightly differs from that derived 
in | Ardavan et al.| |2003 ) because, here, we have eliminated the 
usual step of integrating the first member of equation ( |A20| ) by 
parts to transfer the differentiation with respect to t from the 
delta function onto j. Though making no difference in the sublu- 
minal case, this alternative procedure has the effect, in the super- 
luminal case, of removing a multiplicative factor fi± from 
and instead multiplying the expression for E n by n. 



